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The electric quadrupole moment for the Ad state of **Sr"^, one of the most important 

candidates for an optical clock, has been calculated using the relativistic coupled-cluster theory. The 
result of the calculation is presented and the important many-body contributions are highlighted. 
The calculated electric quadrupole moment is (2.94 ± 0.07)eao, where ao is the Bohr radius and e 
the electronic charge while the measured value is (2.6 ± 0.3)eao- This is so far the most accurate 
determination of the electric quadrupole moment for the above mentioned state. 

PACS numbers: 31.15.Ar, 31.15.Dv,32.30. Jc 



The frequencies at which atoms emit or absorb electro- 
magnetic radiation during a transition can be used for 
defining the basic unit of time 0, S 0- The transi- 
tions that are extremely stable, accurately measurable 
and reproducible can serve as excellent frequency stan- 
dards P, 0] . The current frequency standard is based on 
the ground state hyperfine transition in ^^^Cs which is 
in the microwave regime and has an uncertainty of one 
part in 10^^ 0. However, demands from several areas of 
science and technology have lead to a worldwide search 
for even more accurate clocks in the optical regime. 
The uncertainty of these clocks is expected to be about 
1 part in 10^^ or 10^^ 0. Some of the prominent candi- 
dates that belong to this category are i^^Hg+ j^, ^*Sr+ 
E Jji ^''^Yb"'" [3 etc. Indeed detailed studies on these 
ions will have to be carried out in order to determine 
their suitability for optical frequency standards. In a 
recent article [lOj Gill and Margolis have discussed the 



quadrupole moments of ^^Sr"*", i^^Hg^ and ^^-'^Yb+ have 
been measured in order to determine this shift. However, 
accurate measurements of these shifts are very difficult. 
The calculation of these shifts although very challenging 
could in some cases be more accurate than the measure- 
ments. Such calculations have received relatively less at- 
tention so far. The most rigorous calculation to date has 
been performed by Oskay et.al for i^^Hg"*" using the 
relativistic configuration interaction (RCI) method with 
a multi-configuration Dirac-Fock - extended optimized 
level (MCDF-EOL) orbital basis. We focus on the clock 
transition 5s ^5*1/2 — 4d ^D^n in **Sr+ in the present 



paper. Figure ^ shows the clock transition in ^^Sr+. 
The electric quadrupole moment in the state 4d 
was measured experimentally by Barwood et al. at NPL 
[l5| . Since the ground state 5s does not posses 

any electric quadrupole moment, the contribution to the 
quadrupole shift for the clock frequency comes only from 
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merits of choosing ^**Sr+ as a candidate for an optical the 4c? state. In this letter, we present our rela- 

clock. Till recently, the most accurate measurement of 
an optical frequency was for the clock transition in ^^Sr+ 
which has an uncertainty of 3.4 parts in 10^^ How- 
ever, recently, Oskay et al. ^3 have measured the optical 
frequency of i^^Hg+ to an accuracy of 1.5 parts in 10^^ 
and further improvements are expected [l^ . 

When an atom interacts with an external field, the 
standard frequency may be shifted from the resonant fre- 
quency 01 . The quality of the frequency standard de- 
pends upon the accurate and precise measurement of this 
shift. To minimize or maintain any shift of the clock fre- 
quency, the interaction of the atom with it's surroundings 
must be controlled. Hence, it is important to have a good 
knowledge of these shifts so as to minimize them while 
setting up the frequency standard. Some of these shifts 
are the linear Zeeman shift, quadratic Zeeman shift, 
second-order Stark shift and electric quadrupole shift due 
to the interaction of atomic electric quadrupole moment 
with the gradient of electric field J^] . The largest source 
of uncertainty in frequency shift arises from the electric 
quadrupole shift of the clock transition. The electric 
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FIG. 1: Diagram indicating the clock transition in 
Energy levels are given in cm^^. 



^Sr^ 



tivistic coupled-cluster (RCC) calculation of the electric 
quadrupole moment of ^^Sr+ in the 4(i state. RCC 

is equivalent to all-order relativistic many-body pertur- 
bation theory (RMBPT). 

The details of this theory have been discussed in several 
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articles 0, 0|. Here we shall only give a brief outline. 
Treating the closed shell Dirac-Fock (DF) state |$) as the 
reference state, the exact wave function in RCC theory 
can be expressed as, 



\^) = exp(r) 1$) 



(1) 



where T is the core electron excitation operator. In the 
coupled cluster singles and doubles (CCSD) approxima- 
tion, T can be expressed as the sum of one- and two-body 
excitation operators, i.e. T — T1 + T2, and can be written 
in the second quantized form as 



alalabaJll (2) 



2 /-^ ■^p"-q 
ap abpq 

where and are the amplitudes of the singles and 
double excitation operators respectively. The normal or- 
dered Hamiltonian can be written as, 

HN^H-{^H\<i>)=H-EDF, (3) 
where H is the Dirac-Coulomb Hamiltonian given by 



^ 1 



(4) 



For a single valence system we define the reference state 
as. 



N+l 



)^at|$) 



(5) 



with the particle creation operator a\ . The wave function 
corresponding to this state can be written as. 



N+l 



where, 



)=exp(r){(l + 5„)}K+i) 



(6) 



Sv = + S2v = 2_^ alaygP + alalabaysll , (7) 

v^p bpq 

where S corresponds to the excitation operator in the 
valence sector and stands for valence orbital and sfj 
and s^l are the amplitudes of single and doubles exci- 
tations respectively. Details concerning the evaluation 
of the closed and open shell amplitudes have been dis- 
cussed earlier T^. Triple excitations are included in our 
open shell RCC amplitude calculations in an approximate 
way(CCSD(T))[l3. 

The expectation value of any operator O with respect 
to the state is given by. 



(0) = 



_ {l + ^t}0{l + g}|$A^+l) 

~ ($JV+i| {1 + S-t} exp(Tt) exp(T) {1 + S} 
where O = exp(r1')0 exp(r). 
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FIG. 2: The diagrams (a) and (c) are subsets of dressed pair 
correlation (DPC) diagrams. Diagram (b) is one of the direct 
dressed core-polarization (DCP) diagram. 



The first few terms in the above expectation value 
can be identified as O, OSi, OS2, SlOSi etc; are re- 
ferred to as dressed Dirac-Fock (DDF) , dressed pair cor- 
relation (DPC) (Fig|2Ia)) and dressed core polarization 
(DCP)(Fig|2fb)) respectively. We use the term 'dressed' 
because the operator O includes the effects of the core ex- 
citation operator T. Among the above, we can identify 
few other terms which play crucial role in determining 
the correlation effects. One of those terms is SlOSi+c.c 
((FiglH^c)) which is called as dressed higher order pair 
correlation (DHOPC) since it directly involves the corre- 
lation between a pair of electrons. In table HII individual 
contributions from these diagrams are listed. 

The orbitals used in the present work are expanded in 
terms of a finite basis set comprising of Gaussian type 
orbitals (GTO) 



F^,k{r) = r'=exp(-air2). 



(9) 



with = 0, 1, 2 • • • for s,p, c?, • • • type functions, respec- 
tively. The exponents are determined by the even tem- 
pering condition [gj 



(10) 



The starting point of the computation is the generation 
of the Dirac-Fock (DF) orbitals which are defined on 
a radial grid of the form 



n = ro [exp(?; - l)h - 1] 



(11) 



with the freedom of choosing the parameters tq and h. 
All DF orbitals are generated using a two parameter 
Fermi nuclear distribution 



Po 



1 + exp((r — c)/a) ' 



(12) 



where the parameter c is the half charge radius and a 
is related to skin thickness, defined as the interval of 
the nuclear thickness in which the nuclear charge density 
falls from near one to near zero. The interaction of the 
atomic quadrupole moment with the external electric- 
field gradient is analogous to the interaction of a nuclear 
quadrupole moment with the electric fields generated by 
the atomic electrons inside the nucleus. In the presence 
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of the electric field, this gives rise to an energy shift by 
coupling with the gradient of the electric field. Thus the 
treatment of electric quadrupole moment is analogous 
to the nuclear counterpart. The quadrupole moment 
of an atomic state [^'(7, J, A/)} is defined as the diago- 
nal matrix element of the quadrupole operator with the 
maximum value Mj, given by, 

© = (*(7 J J)|e,,|*(7 J J)). (13) 

Here 7 is an additional quantum number which dis- 
tinguishes the initial and final states. The electric 
quadrupole operator in terms of the electronic coordi- 
nates is given by, 

j 

I 

(J/II = O/ll Cl^^ m ]^Jdr |j2(fcr) -I- 



In Ea. (|15|l . the subscripts / and i correspond to the final 
and initial states respectively and jm{kr) is the Bessel 
function of order to; P and Q are the radial part of the 
large and small components of the single particle Dirac- 
Fock wavefunctions respectively, ji is the total angular 
momentum and Ki is the relativistic angular momentum 
quantum number for the ith electron. The angular factor 
is given by 

O/ll = (-l)0/+i/2)^(2j^ + 1)7(2^ 

><(-V2 l/2)-(^'^'^') (16) 

where 

nd k I') — I ^ ii I + k + I' even 
^ ' ^ ' 1 otherwise 

I and k being the orbital angular momentum and the 
rank respectively. 

Finally using the Wigner Eckart theorem we define the 
electric quadrupole moment in terms of the reduced ma- 
trix elements as 

0/1 e(^) m = (-1)^/-'"/ I ) {jfW 6(2) 

(17) 

In table we have presented the details of the basis 
functions used in this calculation and in table ^ con- 
tributions from different many body terms. The value 
of in the 4d state measured experimentally 

is (2.6 ± 0.3)eao 0, where e is the electronic charge 
and Uq is the Bohr radius. Our calculated value for the 
4d ^£'5/2 stretched state is (2.94±0.07)eag. We have esti- 
mated the error incurred in our present work, by taking 



where the sum is over all the electrons and z is the co- 
ordinate of the jth electron. To calculate the quantity 
we express the quadrupole operator in its single particle 
form as 



e<^) = E'^L^^ (14) 



and the single particle reduced matrix element is ex- 
pressed as |22(1 



+h{krr-^^^ ^ (P«,Q«. +Ok,P«0|. (15) 

I 

the difference between our RCC calculations with sin- 
gles, doubles as well as the most important triple excita- 
tions (CCSD(T)) and only single and double excitations 
(CCSD). A non-rclativistic Hartree-Fock (HF) determi- 
nation resulted in = S.OSeap |l5|. A subsequent cal- 
culation based on RCI with MCDF-EOL orbital basis 
yielded = 3.02eag ^2|- that calculation, correla- 
tion effects arising from a subset of the terms Si , Ti for 
single excitations and 5*2 and T2 for double excitations 
were considered, where and 5*2 are the cluster oper- 
ators representing single excitations from the valence 5s 
orbital to a virtual orbital and double excitations from 
the valence 5s and the core {4s, 4p, 3d} orbitals, with 
atmost one excitation from the core, respectively. In our 
calculation, in addition to these effects, the effects arising 
from the non-linear terms like T|, T1T2, etc. have been 
included. In the framework of CCSD theory, the single 
and double excitations have been treated to all orders in 
electron correlation including excitations from the entire 
core. This amounts to a more rigorous treatment of elec- 
tron correlation in comparison to the previous calculation 
performed using the RCI method.. 

It is clear from the table that the DDF contribu- 
tion is the largest. The leading correlation contribution 
comes from the DPC effects and the DCP effects are an 
order of magnitude smaller. This can be understood from 
the DPC diagram (Figl^a)) which has a valence electron 
in the 4^5/2 state. Hence the dominant contribution to 
the electric quadrupole moment arises from the overlap 
between virtual ^5/2 orbitals and the valence, owing to 
the fact that 6*1 is an operator of rank and the electric 
quadrupole matrix elements for the valence 4^5/2 and 
the diffuse virtual orbitals are substantial. On the 
other hand, in the DCP diagram (Fig|2Ib)), the matrix 
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TABLE I: No. of basis functions used to generate the even 
tempered Dirac-Fock orbitals and the corresponding value of 
ao and (3 used. 



51/2 Pl/2 P3/2 ^3/2 ^5/2 fb/2 jl 12 97/2 ff9/2 

Number of basis 38 35 35 30 30 25 25 20 20" 
Qo(xlO"^) 525 525 525 425 425 427 427 425 425 

P 2.33 2.33 2.33 2.13 2.13 2.13 2.13 1.98 1.98 

Active holes 10 10 10 11 11 8 8 6 6 
Active particles 433110000 



TABLE IT: Contributions of the electric quadrupole mo- 
ment in atomic units from different many-body effects in 
the CCSD(T) calculation. The terms like DDF, DCP, DPC, 
DHOPC are explained in the text. The remaining terms in 
Eq|H|are referred to as 'others'. 



DDF DPC DCP DHOPC Others Total 



3.4963 -0.4306 -0.0642 0.0353 -0.0271 2.94 



element of the same operator could also involve the less 
diffuse s or p orbitals. Hence, for a property like the elec- 
tric quadrupole moment, whose magnitude depends on 



the square of the radial distance from the nucleus, this 
trend seems reasonable, whereas for properties like hy- 
perfine interaction which is sensitive to the near nuclear 
region, the trend is just the opposite for the d states 
[231. As expected, the contribution of the DHOPC ef- 
fect i.e., S\OSi (FiglSfc)) is relatively important as it 
involves an electric quadrupole matrix element between 
the valence 4c?5 /2 and a virtual ^5 /2 orbital. Unlike many 
other properties particularly the hyperfine interactions, 
an all-order determination of this diagram is essential for 
obtaining an accurate value of the electric quadrupole 
moment. One of the strengths of RCC theory is that it 
can evaluate such diagrams to all orders in the residual 
Coulomb interaction. 

In conclusion, we have performed an ah initio calcula- 
tion of the electric quadrupole moment for the 4d 
state of ®*Sr+ to an accuracy of less than 2.5 % using 
the RCC theory. Evaluation of correlation effects to all 
orders as well as the inclusion of the dominant triple ex- 
citations in our calculation was crucial in achieving this 
accuracy. This is currently the most accurate determi- 
nation of electric quadrupole moment for this state in 
Sr^ . Our result will lead to a better quantitative under- 
standing of the electric quadrupole shift of the resonance 
frequency of the clock transition in **Sr+. 
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